Order of Degrees of Freedom in Underdoped Manganites 
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Using the degenerate double exchange Hamiltonian with on-site Coulomb interactions we show 
how the spin, charge, and orbital state degrees of freedom vary in Lai-ajCa^MnOa for dopings in 
between x — 4 and x = 5. With the ordered configuration of the system we investigate the band 
structure for different values of the doping. 



A large family of Mott insulators are transition metal 
oxides with a perovskite structure. Recently, manganites, 
with their particular colossal magnetoresistance, have 
been the subject of intense s tudjii*2* 3 ^'£*£*£'&. Interest- 
ingly, these compounds are very rarely both magnetic 
and ferroelectric, which would allow for technological ap- 
plications. It has been shown that charge order results 
in a lack of inversion symmetry that gives rise to an elec- 
tric dipole momenliii^. In the present paper we extend 
the previous study on the interaction between charge, 
orbital, and spin ordering in underdoped manganitesiiS. 

The main complications in finding the charge, orbital, 
and spin order are the anisotropy of the orbital states 
per site, and the required dependence of the spin cant- 
ing on the doping. The fivefold degeneracy of the d or- 
bital states of the Mn ion is partially lifted by the crystal 
field generated by the surrounding oxygen atoms in the 
perovskite structure. The result is twofold degenerate e g 
orbital states, d x 2_ y 2 and d^ z i- r 2, and threefold degener- 
ate t2 g orbital states, d xy , d yz , and d zx . The degeneracy 
of the e g states is further lifted by the lattice distortions, 
or by the on-site Coulomb interactions. The difference 
is that the lattice distortions lift the degeneracy com- 
pletely by lower the symmetry of the cubic crystal which 
lowers the energy of the d 3z 2_ r 2 state while raising that 
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The on-site Coulomb interactions make the 



system quasidegenerate because the degeneracy is only 
lifted after either of the two orbital states is first occu- 
pied. In an array of manganese ions the orbital states 
in neighboring ions may be rotated with respect to each 
other. After choosing a basis we may characterize the 
orbital states of every ion by an angle 2 . 

The electronic configuration of each ion in the array 
may also vary. The Mn 3+ ion in the parent compound 
of interest, LaMnC>3, has an open shell configuration in 
which three electrons are in the ti g orbitals and one in ei- 
ther of the two e g orbitals. The tig electrons are strongly 
stabilized by the crystal field, and have a lower hybridiza- 
tion with the 2p states of the oxygen atoms, which re- 
sults in the electrons being localized^. When doping 
with a divalent ion that introduces holes, for instance 
Lai-zCa^MnOa, some sites change to a Mn 4+ closed 
shell configuration in which the e g orbitals are empty. 

There are three interactions that reproduce the prop- 
erties of the solids in discussion: the double-exchange, 
the superexchange, and the on-site Coulomb interaction. 
The double-exchange has an energy scale given by the 
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FIG. 1: The lattice with the corresponding unit cell for the 
Zener polaron phase phase <f> = 7r/2 (top), the intermediate 
phase cj> — 7r /4 (middle), and the CE phase = (bottom). 
The rotation of the spins is depicted for each spin canting 
<j>. The Zener polaron phase shows dimers, with blue and red 
dimers AFM coupled within a diagonal. The intermediate 
phase shows the superposition of the Zener polaron and the 
CE phase. The CE phase is composed of AFM coupled zig-zag 
chains 



hopping amplitude t. We consider t = 0.622 eV, which 
is found in the literature^. It leads to the formation of 
bands and to the conduction of electrons. The superex- 
change is a magnetic interaction between localized spins 
of neighboring sites with an energy scale J. The on-site 
Coulomb interaction is the electrostatic interaction be- 
tween electrons in different e g orbital states in the same 
Mn ion. Its energy scale is U. In the Mott insulator the 
otherwise itinerant electrons are localized due to U >> t. 

We consider two limiting cases of Lai_ K Ca K Mn03: the 
Zener polaron phase at x — 0.4 doping, and the CE 
phase at x = 0.5. Both phases are known to be insu- 
lating. As shown in Fig. ^ the unit cell is composed of 
four sites, with one Mn ion each. The CE phase is made 
up of zig-zag chains that have ferromagnetic spin order 
within the chain, and that are antiferromagnetically cou- 
pled with the neighboring chains. This phase presents 
the checkerboard charge order in which sites alternate 
between Mn 3+ and Mn 4+ . The orbital order is d 3z 2_ r 2 
and d x 2_ y 2 in the corner sites (CI and C2), d$ x 2_ r 2 in 
the horizontal bridge (B2), and d3j,2_ r 2 in the vertical 
bridge (BI). The Zener polaron phase is made up of 
dimers with equivalent sites: Bl and CI, and B2 and 
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C2 become equivalent. The spins within a dimer are 
parallel, while spins in between neighboring dimers with 
the same crystallographic direction couple antiferromag- 
netically. Perpendicular dimers have perpendicular spin 
orientations with respect to each other. The charge in 
a dimer is delocalized in between the equivalent sites, 
which results in a bond-centered charge order—. 
The Hamiltonian of the system is given by 

h = t (xilxiK^rf/ 

+ Jaf Y,&'$j + U E «'> (!) 



<y> 



where i,j are site labels, and a, a' denote e g orbital 
states. The sum < ij > is nearest neighbor sites. The 
densities n\ are per orbital per site, and S% are the 
localized (classical) t2 g spins in nearest neighbor sites. 
We consider an antiferromagnetic superexchange inter- 
action, Jaf- The kinetic term contains the fermionic 
creation and annihilation operators, c] and Cj i<7 i , in be- 
tween orbitals of neighboring sites, and the symmetry 
matrix element that takes into account the overlap 
of the orbitals involved and the crystallographic direc- 
tion of the hopping. We consider the overlap of the spin 
wave functions^, (Xi\Xj)> f° r the hopping between near- 
est neighbor sites. It depends on the relative direction 
between the two spins, | (Xi\Xj) I — cos(</>y/2), where faj 
is the relative angle between two spins. The hoping is 
largest for parallel spins, while antiparallel spins cannot 
hop. 

We consider the zero temperature limit with strong 
ferromagnetic Hund's coupling between localized t<i g and 
itinerant e g electrons, Jn/t — > oo. The antiferromag- 
netic interaction becomes Jaf E<y> \S\ 2 cos fcj . With 
respect to the symmetry matrix, we choose our orbital 



basis states to be a 
case we have 
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that takes into account the overlaps of the e g orbitals 
involved when hoping along the y-axisA&. The hopping 
along the cc-axis has a phase with respect to the y-axis 

= -a/3. After 



which results in T 
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choosing a basis we may refer the orbital states of every 
ion, on and ft, to the reference basis with the transfor- 
mation 
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Finally, we take the on-site Coulomb interaction to be 
U/t = 4 and use the mean-field approximation of the 
densities < nf > mf =< cJ )tr Cj )(T >. 

After decoupling the on-site Coulomb interaction, 



first calculate the charge order (CO), orbital order (OO), 
and magnetic order (SO) for different fixed values of the 
doping, varying the spin canting angle from <pij — (CE 
phase) to <£y — ir/2 (Zener polarons). Taking the rela- 
tive angle between spins to change equally we drop the 
site indices and in what follows use <j>. We also bear in 
mind that the magnetic correlation between some spins 
does not change. The system is taken into account by 
rotating the orbital states per site to our working bases, 
a and ft with the transformation 



nf* sin 2 (&/2) + nf cos 2 (&/2) 
nf* cos 2 (&/2)+nf sin 2 (&/2) 



(4) 



before solving numerically. The self consistency equa- 
tions, 
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(5) 



ft 

TV- 



n"nf ~< nf > nf + nf < nf > 



< nf >< xi; >, we 



tan- 1 (2<ct a c i ^>/(n?-nf)) 
nf cos 2 (&/2)+ < c[ a c hl3 > sin(Ci) 
+ nfsin 2 (6/2) 

nf sin 2 (^/2)- < c\ a c %Si > sin(&) 
+ nf cos 2 (^/2), 



enforce that the rotation of the calculated densities, hf, 
from the working basis back to the on-site basis results 
in a diagonal on-site Coulomb interaction. 

Next, we calculate the total energy of the system and 
identify the ground state energy, Eq(<P), for fixed doping. 
The total energy is composed of a term linear in the 
mean- field value of the densities, and a term quadratic 
in the mean-field value of the densities. The linear term, 
namely, the kinetic term, increases monotonously when 
changing the spin canting from one limiting value of the 
doping to another. On the other hand, the quadratic 
term has a clear minimum that moves from one limiting 
doping to the other as the spin canting changes, shown 
in Fig. It is this term that, when added to the linear 
term, gives a minimum to the total energy. Hence, the 
on-site Coulomb interactions stabilize the system. We 
found that for larger ratios of U/t the range of stability 
goes below 0.4 doping, but the validity of the numerical 
calculation is lost as we go beyond the ratio for which 
mean-field is valid. 

Following the determination of Eq(4>) for different dop- 
ing, we stablish the relation between the spin canting <p 
and the doping x by fitting the curve <fi(x) given by the 
energy minima. We use 0(x) to find the ground state 
energy, Eq[<P(x)], as the doping x changes continuously 
from 0.4 to 0.5, i.e., the stable configurations of charge, 
orbit, and spin order. 

With respect to our results for the orbital order, we 
first need to understand the difference between corner 
and bridge orbital states for the CE phase, £ = ±7r/3. We 
bear in mind^ that pairs of orthogonal orbital states differ 
by an angle of n, |£ > and |£ + n >. Since our basis or- 
bital states are chosen to be |0 >=3z -r and \tt >=x 2 -y 2 , 
and the orbital states at the bridge sites having nonzero 
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FIG. 2: Contribution to the energy from the quadratic term 
in the mean-field value of the densities for fixed doping as 
a function of spin canting for U/t = 4. For different doping 
values between the Zener polaron phase (top) to the CE phase 
(bottom) with Aa; = 0.02 increments. Notice the change of 
the minimum value for each doping as a function of (j>. 



occupation are \ir/3 >=3y 2 -r 2 and | — tt/3 >=3x 2 -r 2 , we 
conclude that the reference orbital state for £ = ±7r/3 is 
\ir >=x 2 -y 2 . This allows us, on the one hand, to interpret 
the found mean field densities of orbital states per site 
by noticing that the labels a and (3 at the bridge sites 
are interchanged with respect to the corner sites. That is 
to say, we have to take into account the phase difference 
in orbital states when the charge moves from a |£ + it > 
orbital state to a |£ > one. On the other hand, we can 
expand on the idea that the on-site Coulomb interaction 
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FIG. 3: Charge disproportion Aq and spin canting <f> as a func- 
tion of doping. The CE phase shows the largest disproportion 
between sites, corresponding to the checkerboard charge or- 
der, while the Zener polaron has Aq = corresponding to 
charge ordered on the bond between bridge and corner sites. 
For the spin canting we see a continuous rotation between CE 
and Zener polaron phase. 



pushes the charge from the corner into the bridge sites™ . 
Previously, no reference to the specific orbitals involved 
in the hopping was made. We conclude that the hop- 
ping takes place from the x 2 -y 2 orbital states at the cor- 
ner sites unto the 3y 2 -r 2 and 3x 2 -r 2 sites at the bridges 
given the required rotation of orbital states. Since the 
difference is exactly |vr/3| it seems to suggest that all 
the hopping between bridge and corner sites take place 
in the above mentioned orbitals, and none from the 3z 2 - 
r 2 orbital state at the corners into the bridges. This is 
in agreement with the findings of other computational 
schemes^. 

Interestingly, unlike the change of density per orbital 
state, the angle £j does not change linearly with the dop- 
ing for both bridge and corner sites orbital states. As 
seen in Fig. the orbital states at the bridge sites re- 
main almost constant between x = 0.45 and x = 0.5. For 
the Zener polaron phase we find that the orbital states 
per site for each dimer is the same, as required for the 
charge to be delocalized and centered on the bondi. In 
this limit we find that the orbital states are defined by 
£ = ±7r/4, with the angular difference between the or- 
bital state in dimers along the y-axis with respect to those 
in the cc-axis equal to \n/2\. Thus, the orbital states be- 
tween horizontal and vertical dimers are orthogonal with 
respect to each other. The result is two sublattices of 
dimers in diagonal stripes that have orthogonal spin and 
orbital order with respect to each other. 

Concerning the charge order, our results reproduce the 
checkerboard charge order of the CE phase, and the pre- 
dicted charge order for the Zener polaron phase, which 
is on the bond±. This can be seen in Fig. [21 where the 
charge disproportionation between two neighboring sites 
is zero for the Zener polaron phase, and maximum for 
the CE phase. Previously, a similar curve had been pro- 
posed for U/t = 1, were the largest valence difference 
was less than Aq = O.Ie at the CE phasei. We find 
for U/t — 4 a charge disproportionation slightly greater 
than Aq = 0.2e. As already known, a stronger on-site 
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FIG. 4: The charge order (top) and orbital order (bottom) 
for different doping. 
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FIG. 5: Bands between points of high symmetry of the re- 
ciprocal space lattice for U/t — 4. Notice the changes of 
the energy gap between valence and conduction bands, and 
the closely packed and almost dispersionless bands above the 
conduction bands. 



Coulomb interaction pushes more charge from the cor- 
related corner sites unto the uncorrelated bridge sites 
of the CE phase. This result is far from the one us- 
ing exact diagonalization but well within the mean-field 
approximation^. We also present in this graph our re- 
sults for the dependence of the spin canting <f> on the 
doping. 

The density of particles per orbital state per site is 
plotted in Fig. At x = 0.5 the corner sites have 
more charge in the orbital state d x 2_ y 2 (C 13 ), than in the 
d 3z 2_ r 2 (C Q ). On the bridge sites all the charge is in the 
orbital states d3 a .2_ T .a/d3, ( 2_ r a (B 13 ), while the B a orbital 
is empty. Hence, the bridge sites are uncorrelated, as al- 
ready showniS. Notice how the aforementioned exchange 
of labels for the orbital states is crucial in interpreting 
the densities properly. The densities change linearly as 
a function of doping. As the doping decreases from the 
CE phase we see two main effects: the densities C° and 
B 13 become equal as well as C 13 and B a , and the density 
B' 3 becomes nonzero. Hence, the corner and bridge sites 
at x = 0.4 become equivalent, resulting in Zener polaron 
dimers. The equivalence between sites results in a der- 
ealization of the charge in between the two, that leads 
to the charge being centered on the bond on average, as 
previously mentioned. 

After finding the stable configuration of the system 
we calculate the band structure for different spin cant- 
ing values, shown in Fig. [S] Doing basic band structure 
analysis^ we see for the CE phase that the dispersionless 
bands along T — (tt, 0) and (tt, tt) — (0, tt) in the reciprocal 
lattice indeed show antiferromagnetically coupled zig-zag 
chains. Moreover, the bands are degenerate along this di- 
rection. We also see how the energy gap at the T point be- 
tween the valence bands and the empty conduction bands 
is the largest with respect to the other spin canting bands 
shown. For half-filling all the valence bands are filled and 



the system is an insulator. In addition, we also notice 
that four of the bands are above the conduction bands: 
the strong correlation between on-site orbitals for U ^ 
makes four bands energetically inaccessible by shifting 
them above the conduction bands. This makes the sys- 
tem effectively monovalent. In the limit U/t — > these 
bands become fourfold degenerate zero energy bands. It 
was previously argued that while in the CE phase one 
could dope these nondispersive bands causing no change 
in energjiiS, thus possibly occupying both orbital states 
per site. We add to it that, as seen in Fig. [f)| this is 
clearly not the case as soon as (f> ^ 0. For spin cant- 
ing (j) = 7r/4, we see that the bands along T — (it, 0) and 
(tt, tt) — (0, tt) are now active, and cross the Fermi energy 
^0_7r/4y^ _ _q g4g a i on g t ne segment r — (tt, tt). Further, 
their degeneracy has been lifted by the opening of a small 
energy gap at the boundary of the Brillouin zone. On the 
other hand, the energy gap at the T point between the 
valence and conduction bands closes. The energy gap 
splitting occupied and empty bands in the CE phase in 
the limiting case U — ► was found^ to be Ay^ /t = 1. 
For U j= the gap reduces to A CE /t = 0.83, which 
is A CE — 0.51 eV. That the energy gap may be small 
enough to allow for double-exchange has been discussed 
experimentally 8 . Finally, for the Zener polaron phase 
at <t> = 7r /2 we see that the system is practically dis- 
persionless along the boundaries of the reciprocal lattice 
(tt, 0) — (tt, tt). Also interesting is that the bands now 
cross each other along the path T — (tt, tt), and cross the 

fermi energy e^ - ^^ 2 /t = —0.835. The energy gap further 
closes to A zp /t = 0.57, which is A zp = 0.35 eV. 

Finally, the density of states, p(e), is calculated. It has 
sharp narrow peaks corresponding to the bands above 
the conduction bands. Focusing on the density of states 
of the valence electrons, Fig. [7| we notice that the Fermi 
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FIG. 6: Total energy for the limiting case U — » 0. Notice how 
the possibility to fill zero energy states changes as soon as the 
spin canting becomes nonzero. 
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FIG. 7: Density of states corresponding to the valence bands. 
The dotted lines correspond to the values of the Fermi energy 
of the insulating CE phase (right) , intermediate spin canting 
4> = 7r/4 (middle), and Zener polaron phase (left). 



energy of the CE phase is delocalized at the middle of the 
gap between valence and conduction states, giving rise 
to an insulating state. Interestingly, for the intermediate 
phase at cj> = 7r/4 the Fermi energy coincides with a van 
Hove singularity. This is due to the point (tt, 0) of the 
reciprocal lattice, where the group velocity vanishes. 

We conclude that in undcrdoped Lai-aCa^MnOa for 
doping in between x = 0.4 and x — 0.5 the effect of 
the on-site Coulomb interactions is to make the solid ef- 
fectively monovalent, and to stabilize the orbital, charge 
and spin configuration. We have also shown how the 
orbital, charge and spin order change when going from 
the Zener polaron phase to the CE phase. Regarding 
the band structure, it was proven that the energy gap 
closes as the spin canting changes from zero to 7r/2, and 
is smaller than the hopping amplitude, possibly allow- 
ing for double-exchange processes. It has been proposed 
elsewhere^ that the inclusion of the oxygens in the calcu- 
lations results in the CE phase being made up of oxygen 
stripes. Since we are considering atomic orbital wave 
functions our description of charge order in the bond is 
an assumption that does not allow us to compare our re- 
sult with this theoretical work. Our results show that the 
magnetic correlation between spins within a dimer is al- 
ways ferromagnetic and does not change for the dopings 
being discussed. We suggest that within the limitations 
of our model both results are related. We do not find 
any agreement with the results in which a different con- 
figuration of the CE phase was proposed^. Future work 
should consider thermal effects. 
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